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'. Abstract 

QQ ■ As a step towards understanding multi-quark systems abundant in nature we 

CN . construct a model that reproduces the binding energies of static four-quark 

. systems. These energies have been calculated using SU(2) lattice gauge theory 

J> I for a set of six different geometries representative of the general case. The 



(N 



model is based on ground and excited state two-body potentials and 
multi-quark interaction terms. 

o> ■ 

O : PACS numbers: 11.15.Ha, 12.38.Gc, 13.75.-n, 24.85,+p 

^ ; 1 Introduction 

^ ■ The significant progress in lattice QCD has so far been restricted to systems with 

only a few quarks - up to three in most cases. However, in particle and nuclear 
physics there is also considerable interest in few- and multi- /iadron systems be- 
^ I ginning with the possibility of bound KK states. This prompts one to ask what 

I lattice QCD can say about these more complex quark systems. It seems unlikely 

that, in the foreseeable future, the techniques of lattice QCD can be developed 
sufficiently to tackle these problems directly. Therefore, we are reduced to con- 
structing models that can explain the lattice results for the simplest multi-quark 
systems in a way that can be readily extended to more complicated cases. With 
this in mind, the Helsinki group in Refs. ||l], ||, |^ has calculated the energies 
of four quarks on a lattice in various geometries - namely, four quarks at the 
corners of rectangles, tetrahedra and other geometries as shown in Fig. |I[ This 
is taken to be our 'experimental' data, which is to be explained by some model. 
Hopefully, this set of geometries is general enough for the model to also explain 
the energies of geometries in between those actually considered. Our philosophy 
is that, if any geometry cannot be fitted, then the model fails, since then there 
is no reason to expect configurations not checked explicitly to be fitted. 
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^E-mail: Petrus@hip.fi 
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Due to limitations of computing resources the 'experimental' data is not for full 
QCD. Instead of SU(3) we use SU(2), which saves about an order of magnitude 
in computer time. The less than 10% quantitative differences observed in the 
glueball spectrum and topological susceptibility of these theories suggest that the 
relevant qualitative features are preserved [^]. Secondly, since the quarks are fixed 
in some geometry, they are static. Therefore this discussion is most relevant to the 
heaviest quarks, in analogy to the static potentials used successfully to describe 
quarkonia. The static approximation is currently being partially removed by 
applying the techniques of Ref. [^] to a system of two B mesons with only the b 
quarks static. A third approximation is that no quark-pair creation is allowed the 
so-called quenched approximation. The effects from this approximation have been 
found to be of the same magnitude as the use of SU(2) instead of SU(3), apart 
from string breaking effects at distances larger than we simulate. We concentrate 
here on attempting to understand the best simulation data available, i.e. for a 
system with four static SU(2) quarks in the quenched approximation. 

Any model that can be extended to multi-quark systems must presumably treat 
only the quark degrees of freedom explicitly - with the gluon degrees of freedom 
entering only implicitly. This is the same philosophy as used, with much success, 
for interacting multi-nucleon systems. In the latter, the meson fields generating 
the interactions lead to effective internucleon potentials, which take the form of 
two-nucleon potentials and, to a lesser extent, three- and four-nucleon potentials. 
However, it must be remembered that multi-quark and multi-nucleon systems are 
very different, since in the former the underlying gluon fields are strongly self- 
interacting. It is, therefore, not at all clear that any effective model defined in 
terms of only quark degrees of freedom will be successful. It is the purpose of this 
article to see to what extent such a model can be developed for the four-quark 
case. If this trial fails, then there is no point in expecting it to work in even more 
complex quark problems i.e. a successful model for four quarks is necessary but 
not sufficient before considering any extension to even more quarks. 

In Section 2 the model is introduced. The results are given in Section 3 and 
concluding remarks presented in Section 4. 

2 The model 

The basis A,B,... of the model contains pairings of the quarks in all possible 
ways. A normalization matrix N contains the overlaps of any two states, Naa =< 
A\A >, Nab =< A\B >,..., while a potential matrix V for an interaction V 
is defined as Vaa =< A\V\A >, Vab =< A\V\B >,.... The object is then to 
compare the eigenvalues 

[V - E{A)N] = (1) 
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with the lattice results - the success or failure of the model being to what extent 
the two can be made to agree. 

As we only consider SU(2), there is no distinction between the group properties 
of quarks (g) and antiquarks (g). As shown in Fig. |l|, four such quarks can then 
be paired in three different ways 

^ = (?i 93) (92^4), B = {qiq2){,q?.qA) and C = (qm) iq2q3) , (2) 

where each (qiqj) is a colour singlet. These three basis states are not orthogonal 
to each other and have in the weak coupling limit the condition 

\A + B + C >=0. (3) 

Since < A\A >=< B\B >=< C\C >= 1, we get in this limit the equalities 
< A\B >=< B\C >=< A\C >= -1/2. 

The form of the potential is motivated by perturbative one-gluon exchange i.e. 

i<j 

where, for the ground state, we use the Cornell form of parameterization of the 
static two-quark potential 

Vij = -e/vij + hsTij + c (5) 
fitted to measured lattice values. Then Vaa = "^13 + "^24 etc. and 

Vab = Vba =< A\V\B >= -i (i;i3 + V24 + vu + V34 - vu - V23) ■ (6) 

Here it should be noted that using this perturbative form of the potential without 
a multi-quark interaction term is unacceptable because of the resulting unphysical 
long range Van der Waals forces. When moving to stronger couplings the overlap 
between any two states should decrease, becoming zero in the strong coupling (or 
large distance) limit. We take this effect into account by introducing a factor /, 
defined as < A\B >= — //2, which decreases with increasing separation of the 
quarks. 

For internal consistency, this same factor must also multiply Vab in Eq- i ^ 
otherwise the eigenvalues would depend on the self-energy term c in Eq. ^. A 
reasonable parameterization for the multi-quark interaction term / is 

/ = exp{-bskfS), (7) 

where bg is the string energy of Eq. ^, 5* the minimum area bounded by the 
straight lines connecting the quarks and kf a parameter to be fitted. 
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This form of the parameterization was originally motivated by strong coupling 
ideas 0, If. More complicated parameterizations, such as / = foexp^—bgkfS — 
\/b^kpP) with P the perimeter of S and fo,kp new constants to be fitted, have 
been suggested. However, the new parameters were found to measure lattice 
artefacts in a continuum extrapolation, where /o — ^ 1, fcp — > |0. The deter- 
mination of a minimal area 5* is a laborious task [|10] and thus impractical when 



many such areas are needed. Therefore, we use a formula that is readily evalu- 
ated for any geometry; the area is simply taken to be the average of the sum of 
the four triangular areas defined by the positions of the four quarks i.e. the faces 
of the tetrahedron. For example, in the notation of Eq. ^, the appropriate area 
S{AB) for / is 

S{AB) = 0.5[S(431) + S(432) + S(123) + 5(124)], (8) 

where S{ijk) is the area of the triangle with corners at i,j and k. For planar 
geometries this simply reduces to the expected area, while for non-planar cases 
this is only an approximation to S{AB). 

Of the six geometries considered here the tetrahedron presents the most challeng- 
ing data. Other geometries could be fitted qualitatively with the above model 
using only the one free parameter kf in Eq. 0. However, the regular tetrahedron 
has maximal degeneracy as all quarks are at an equal distance from each other. 
Its energies have the interesting feature that the lowest state is doubly degen- 
erate and becomes more bound as the tetrahedron increases in size, as noted in 



Ref. |TT|. This is opposite to what happens with squares, where the magnitude 
of £^0 decreases with increasing size of the system. This indicates that there is 
coupling to some higher state (s) that becomes more effective as the size increases 
and suggests that these higher states contain gluon excitation with respect to 
the basic A,B,C configurations. Therefore, as proposed in Ref. 0, we include 
states, denoted with a prime, which describe the same quark partitions but with 
the gluonic field in an excited state. These excited states could be excitations 
of the two-body potential or flux conflgurations where the four quarks form a 
singlet. In the latter case the formalism below needs to be slightly modifled. If 
the excited state is taken as the lowest two-body excitation, for which the gluonic 
field has the symmetry of the Eu representation of the lattice symmetry group 
D^h, we have 

A* = {(li(l?i)EM2qi)Eu etc. (9) 

Because the Eu state is an odd parity excitation, A*,B*,C* must contain two 
such states in order to have the same parity as A, B, C . The excitation energy 
of an Eu state over its ground state(y4ig) counterpart is ~ vr/i? for two quarks 
a distance R apart. As said above, when R increases this excitation energy 
decreases making the effect of the A*, B* , C* states more important. 

In analogy to Eq. ^ an antisymmetry condition A* + B* + C* = Q holds for the 
excited states. As the gluonic excitations are orthogonal to the the ground state 
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we also have < A\A* >=< B\B* >=< C\C* >= 0. There are now two more 
functions /"''^ defined as 

< A*\B* >=< A*\C* >=< B*\C* >= -r/2 and 

< A*\B >=< A*\C >= . . . etc ... = -/72. (10) 

Here it is assumed that /"'^ are both dependent on 5* as defined in Eq. p. Since 
involves only the excited states, it is reasonable to expect it has a form similar 
to / in Eq. ^ i.e. 

r = expi-b,KS). (11) 

In the weak coupling limit, from the A* + B* + C* = condition, we expect 

< A\B* >=< B\C* >= = at small distances. This also happens in the 

Isgur-Paton model, as can be seen in Fig. 6 of Ref. p. To take this into account 
we parameterize /" as 

r = {fS + rsS)exp{-bsKS). (12) 

Justification of this is given by the Isgur-Paton calculation and an analysis of the 
correlations between the paths involved in the simulation, which gives the term 
linear in S. When all three parameters are varied in a Minuit fit it is found that 
/o = 0.002(30) is consistent with zero as expected. Therefore, from now on we 
set this parameter to zero. 

For the potential matrix V(f) the diagonal matrix elements, before the lowest 
energy amongst the basis states - usually Vaa - is removed, are 

< A*\V\A* >=v*{13) +v*{2A), etc, 

where v*{ij) is the potential of the Eu state - a quantity also measured on the 
lattice along with the four-quark energies. However, to allow more freedom in 
the following fits, we introduce a parameter bo in 

< A*\V\A* >= Vaa + boV^A, VXa = v*{l?>) + v*{2A) - v{l?>) - v{2A) (13) 

for state A* and analogously for B*,C*. 

If Eu is the most important excitation in the parameter fit, bo should turn out 
to be of the order unity. However, in practice we get 2.3(5), suggesting that 
excitations with higher energies are, perhaps, more relevant. 

The two-quark potentials v{ij) are taken to be more elaborate than the three 
term form of Eq. ^ They are fitted to the lattice data using 



, , 0.255 0.045 , , 

v{rij) = 0.562 + 0.0696rij —. (14) 



5 



Similarly, the excitation of the Eu state is fitted by 



v*{ij) - v{ij) = — 



TT 4.24 3.983 



+ 



(15) 



The extra terms containing and r~j^ are purely for numerical reasons and 
ensure that the fitted values of v{ij) and v*{ij) are, on average, well within 1% 
of the lattice values for all r,. > 2. 



There are two types of off-diagonal element 

< A*\V\B* > and < A\V\B* > . 



(16) 



These require further approximations. Using the Isgur-Paton model with N=l 
as a guide we get qualitatively 



< A*\V\B* >= - 



r 



Vaa + Vbb - Vcc + co 



{VXa + V^b) 



and analogously for the elements < yl*|y|C* >, < B*\V\C* >. 
Likewise, 



:< B*\V\A* > 

(17) 



< A\V\B* >= -4r 
' ' 2 



Vaa + Vbb — Vcc + 



V 



BB 



< A*\V\B >-- 



r 

2 



Vaa + Vbb - Vcc + «o- 



aa 



etc, 



(18) 



where a^, Cq are free parameters, which should have values of order unity if is 
the most relevant excitation. 



In the special case of regular tetrahedra, V reduces to the form 



V = 



Vaa 

-fVAA/2 
-fVAA/2 


-fVAA/2 

Vaa 

-fVAA/2 


-fVAA/2 
-fVAA/2 

Vaa 




-rvj2 
-rvj2 


-rvj2 



-rvj2 


-rvj2 - 
-rvj2 






-rvj2 
-rva/2 


-rvj2 



-rva/2 


-rva/2 
-rva/2 




V 

-rvj2 
-rve/2 


-rvj2 

Vb 

-rvj2 


-rvj2 
-rvj2 



(19) 

where K = Vaa + aoKI^/S , H = Vaa + hoVX^ , K = Vaa + cqVXa- As with all 
geometries 



" 1 

-//2 


-//2 
1 

-//2 


-//2 
-//2 
1 




-/V2 
-/V2 


-/V2 


-/V2 


-/V2 " 
-/V2 





-/72 

_ -/72 


-/V2 


-/V2 


-/V2 

-/72 




1 

-/72 
-/72 


-/72 
1 

-/72 


-/72 

-/72 
1 



(20) 
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The full 6x6 matrix [V — E'N] now breaks into three 2x2 matrices, two of 
which are identical - giving the observed degeneracy. These have the form 



-E(l + //2) -f-(E-Va)/2 

riE - K)/2 -^(1 + r/2) + H + rvj2 



[V - ^N] = 

whereas the third 2x2 matrix giving nondegenerate energies is 



(21) 



[V - ^N] 



-Eil-f) 



nE-Va) 
-E{i - n + H - rv. 



(22) 



3 Results 



In Refs. [H, ^, |TT1 four quark energies have been extracted for a variety of 
geometries using a 16'^ x 32 lattice with (3 = 2.4. This (3 value corresponds to a 
lattice spacing a = 0.119(1) fm. From these energies, one hundred - distributed 
over all measured geometries - are selected for fitting. Configurations containing 
flux links of less than two lattice units were not included because of the strong 
lattice artefacts they contain. The six geometries are shown in Fig. 1. Specif- 
ically, we use 15 Tetrahedra (T), 6 Squares (S), 12 Rectangles (R) (including 
Tilted Rectangles (TR)), 4 Quadrilaterals (Q), 9 Non-Planar (NP) and 4 Linear 
(L). Only the lowest two energies (-Eo,i) from the lattice simulation are used. In 
most three basis simulation had been performed, so that a third energy 

{E2) was in fact available. However, as this state is the highest calculated, it 
is not expected to be very reliable due to the higher excitations it contains. Its 
main purpose was to improve the estimate on Ei by reducing its excited state 
contamination. One might question the connection of these energies with the 
continuum values. This has been answered in Ref. |^, where it was found that 
the binding energies for equal physical sizes essentially stay constant when the 
lattice spacing is made smaller. 

Before commencing a fit, the size of the errors on the above data must be decided. 
The lattice simulation, through the boot-strap method, does indeed produce er- 
rors - statistical ones. However, some estimate must also be added for systematic 
errors. How this is done is somewhat subjective. Here the prescription is to as- 
sume all errors must be at least 0.005 and, also, at least 10% of the eigenvalue 
itself. The former corresponds to about 10%, 1% for the largest values of Eq, Ei 
respectively. 

The above 100 pieces of data were fitted with Minuit - the Migrad option - using 
the seven parameters: kf in Eq. |^, /|, ka in Eq. |12|, kc in Eq. 11, 60 in E<1- 13 and 



ao,co in Eqs. |T8|. The outcome yielded a x^/d.o.f.=1.08 with the values of the 



parameters being given in Table |l[ 
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Parameter 






fs 


kc 


ao 


bo 


Co 


Value 


1.25(6) 


0.54(11) 


0.046(3) 


0.04(20) 


4.4(3) 


2.2(6) 


8.0(4) 



Table 1: The values of the parameters defining the interaction 



Of these values: 

a) The observation that /cc ~ implies that /'^ ~ 1 i.e. the excited configurations 
interact amongst themselves in the way expected from perturbation theory. 

b) The values of bo (along with ao, cq) are somewhat larger than the naively 
expected value of unity. This suggests that higher order effects are important. 

Table |^ shows the contributions to the total x^/ d.o.f. from each of the 12 types of 
data - i.e. from Eo, Ei for the six geometries (T, S", i?, Q, NP, L). For comparison 
the 2x2 basis model {A, B) with / = 1 is also shown. This corresponds to the 
absence of multi-quark interaction. In spite of the frequent use in the literature 
of such models based only on two-body potentials this seems to be a very poor 
choice. 



Basis 


6x6 


2x2 


Geometry 


x\Eo) 




x\Eo) 




T 


0.18 


0.20 


2.4 


21.4 


S 


0.01 


0.11 


6.2 


40.0 


R 


0.17 


0.09 


9.6 


63.9 


Q 


0.09 


0.03 


0.03 


0.25 


NP 


0.08 


0.01 


0.06 


1.53 


L 


0.01 


0.11 


0.01 


0.10 



Table 2: The contributions of the two states {Eo, Ei) of each of the six geometries 
to the total xVd-o.f. of 1.08 (6x6 basis) and 146 (2x2). 



4 Discussion 

The above model with 6 basis states fits well all simulated geometries and confirms 
our earlier work that a 2 x 2 basis model with only two-quark interactions (i.e. 
/ = 1) is not able to even qualitatively account for the data. Areas of further 
study include the actual nature of the higher excitations, which seem to be playing 
an important role. A natural extension of the present static quark model is the 
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application of the model in a dynamic case using the Schrodinger equation. Work 
on simulating such a system is in progress. 
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